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In this article we review the astrophysical application of gravitational microlensing. Af¬ 
ter introducing the history of gravitational lensing, we present the key equations and 
concept of microlensing. The most frequent microlensing events are single-lens events 
and historically it has been used for searching dark matter in the form of compact astro- 
physical halo objects in the Galactic halo. We discuss about the degeneracy problem in 
the parameters of lens and perturbation effects that can partially break the degeneracy 
between the lens parameters. The rest of paper is about the astrophysical applications of 
microlensing. One of the important applications is in the stellar physics by probing the 
surface of source stars in the high magnification microlensing events. The astrometric and 
polarimetric observations will be complimentary for probing the atmosphere and stellar 
spots on the surface of source stars. Finally we discuss about the future projects as space 
based telescopes for parallax and astrometry observations of microlensing events. With 
this project, we would expect to produce a complete stellar and remnant mass function 
and study the structure of Galaxy in term of distribution of stars along our line of sight 
towards the centre of galaxy. 
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1. Introduction 

In this article we will review the astrophysical applications of a specific type of 
gravitational lensing, so-called gravitational microlensing, in the Milky Way stars. 
The bending of light in classical mechanics, when a light ray passes by a massive 
object assuming a mass for the photons is foreseen. The extension of this idea to 
the curved space-time in general relativity leads us to a correct value for the light 
bending which is compatible with the observations. 

While in the gravitational lensing more than one image is produced, in microlens¬ 
ing due to small separation of images, only the magnification effect of a transit lens 
over a background star can be measured. The application of microlensing for explor¬ 
ing the dark matter of Galactic halo answered to one of main astrophysical questions 
about the nature of dark matter in the halo of galaxy. After one decade monitoring 
of Large and Small Magellanic Clouds by three main observational groups, they 
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conclude that dark matter in the halo of Milky Way is not made of Massive Astro- 
physical Compact Halo Objects (MACHOs). 

The other application of microlensing is the possibility of detecting extrasolar 
planets. A planet or a group of planets that orbit around the parent star can play 
the role of multiple lensing in the microlensing events. In this case, parent star with 
the planets around it behave as double or multiple lenses and the result is producing 
caustic lines on the source plane. Once a source star crosses these singular lines, 
a small perturbation on the light curve appears and result is short blip or dip in 
the light curve. Analyzing the light curve reveals characteristics such as distance of 
planet from the parent star and mass of an exoplanet Gravitational microlensing 
compare to the other methods of planet observations, is a unique tool to detect 
planets beyond the snow line as well as planets orbiting in the habitable zone of the 
parent stars.^^®^ The detection of exoplanets by microlensing is a well established 
method and next generation of microlensing experiments with high cadence and 
better photometric precision will detect a larger number of exoplanets.^ 

The new area of gravitational microlensing experiments would be the astrophys- 
ical application of this phenomenon that we will review the details of this subject in 
this paper. The gravitational microlensing is a natural telescope that can probe the 
surface of source stars, located few thousands light years far from us.^ The surface 
of stars are not uniform and due to scattering of light from the atmosphere, for an 
observer, it seems darker at the edge of star and brighter at the central part of it. 
Moreover, stellar spots are defects on the surface of star by means of cooler areas 
compare to the background of star. In the case of high magnification microlensing 
events, lens probes surface of source star and it can encode the surface brightness 
of star in the microlensing light curve. 

The other possibility in the microlensing experiment is the astrometry of images 
as well as polarimetry observations. Both methods will provide extra information 
beside the photometric observations. In the case of astrometry which can be carried 
out by the space-based telescopes or high resolution ground based telescopes, we can 
measure the time variation of centre of light for the source stars during microlensing. 
Astrometric observations directly provides the size of Einstein angle. On the other 
hand polarimetric observations can identify the local polarization on the surface of 
source star and also identify the spots on the stars. This observation will open a new 
window to the stellar studies of remote stars and investigate the physics of stars 
with more details. The other advantage of a space based telescope orbiting around 
Sun is that we can make parallax measurements and with the extra information 
from the astrometry observations, it is possible to identify the mass and distance 
of lenses from the earth. This would be a great possibility for measuring the mass 
function of visible and stellar remnants, as well as the distribution of stars in the 
Milky Way. 

In section ([2]) we will provide a brief history of gravitational lensing and mi¬ 
crolensing. In section o, formalism of gravitational lensing, including the lens 
equation is discussed. Section O introduces calculation for the statistics of mi- 
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crolensing events in different directions of Milky-Way galaxy. In section ([5]), we 
will discuss about the perturbation effects as the parallax and finite-size effects on 
the light curve of microlensing events and applications of these effects for breaking 
the degeneracy between the lens parameters. Also small perturbations on the light 
curve by the effects of stellar spots and compact objects in binary systems will be 
discussed in this section. In section ([61), we review the astrometry and polarimetry 
observations of microlensing events and our prospects for future observations and 
astrophysical applications of polarimetry. The summary is given in section 0. 


2. A brief history of Gravitational Lensing 


In Newtonian gravity, assuming photons are made of tiny particles and moving with 
the speed of light, one can predict the bending of light bundle while passing beside 
a source of gravity. We will use the term of gravitational lens for the gravitational 
deflector of light. For a light ray travelling from a distant source to the observer, 
let us take the minimum distance of the light ray to a point mass lens to be 
Then, the gravitational acceleration perpendicular to the propagation of light ray 
is given by 

GM C 


where x = ct and the time is set to start from — oo and goes to -l-oo. We can 
integrate acceleration in the j/-direction and calculate the velocity of a test particle 
in this direction. The bending angle at each point along the trajectory of light is 
simply giving hy a = Vy/vx and the overall bending angle is 


a = 


2GM 


( 1 ) 


In general relativity which is the extension of special relativity to the accelerating 
frames, the light bending is calculated by the deviation of trajectory of the light rays 
as a null geodesics in the space-time. The standard approach is using Schwarzschild 
metric 

ds^ = —di^{l — -I- r^dO^ + sin^ 9d4i'^, (2) 

r r 

which describes the metric of space-time around a point mass object. We can inves¬ 
tigate the bending of light rays through Fermat principle in which light rays chooses 
the extremum path for traveling from one point to another point. In the case of us¬ 
ing general relativity, the bending angle is twice of that in the Newtonian gravity, 
given in equation O- The prediction for bending of light in General Relativity has 
been confirm by Eddington in 1919. A full story of expedition for 1919 eclipse can 
be found in Ref.^ 

The observability of gravitational lensing where one star plays the role of lensing 
of another background star has been motivated to Einstein by a Czech engineer, 
Rudi W. Mandl in 1936 and Einstein wrote a brief paper by his request in Science.^ 
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However, Einstein noted in his paper ” that is unlike to detect this phenomenon”. 
One year later a Swiss astronomer Fritz Zwicky noted that while the probability 
of observation of lensing of a star by another star is low, however in cosmological 
scales it is likely that galaxies can lens the light of other background galaxies.^ 

In 1986 Bohdan Paczynski returned to the question by Einstein in 1936 whether 
the gravitational lensing of a star by another star inside our galaxy is observable!^ 
He calculate the observability of this effect and proposed this method for indirect 
observation of Massive Astrophysical Compact Halo Objects (MACHOs) in the 
Galactic halo. Assuming that 100% of halo is made of MACHOs, the probability of 
detecting gravitational lensing by monitoring stars in the Large and Small Magel¬ 
lanic Clouds is 10“^. The leading groups of ERO^and MACHO started monitoring 
Magellanic Clouds as well as the centre of Galaxy for searching microlensing events 
and the first microlensing candidate is found in 1993 at the direction of Large Mag¬ 
ellanic Clouds.l^l^ After a decade of monitoring Magellanic Clouds, these groups 
compared the expected number of microlensing events if 100% of halo is made of 
MACHOs, with the observational results. By counting the number of observed mi¬ 
crolensing events, the observation groups came to the conclusion of putting an upper 
limit of 20% contribution of MACHOs in the Galactic halo.^^ 

While the detection of MAGHOs was one of the main goals of gravitational 
microlensing experiments, the application of this method for exploring exoplanets 
and studying the atmosphere of distant stars opened a new window to the astro- 
physical applications of this new tool. The idea of exploring binary lenses by the 
gravitational lensing has been proposed by Moa and Paczynski in 1991. In the later 
years this method become an important tool for detecting planets orbiting stars 
few kilo parsec far from us.l^ They not only could detected earth like planets, also 
their distances from the parent stars were about few astronomical unit, possibility 
of being in the habitable zone. Moreover, due to control on the defection efficiency, 
this method is able to provide correct statistics from the planet census. 

The new surveys as GAIA and space based microlensing telescope^^ as well as 
ground based telescop^^ for the astrometric observation will enable us in near future 
to break the degeneracy between the lens parameters and identify the parameters of 
lenses. One of the main achievement will be the full determination of mass function 
of all compact objects from the brown dwarf to the blackholes. The polarimetry 
observations also will reveal the atmospheric structure of a source star as well as 
the stellar spots of stars. 


3. Formalism of Gravitational Microlensing 

The modern area of microlensing started with the paper of Paczynski in 1986 who 
proposed the observation of lensing of a source star by another star in the Milky Way 
galaxy. In this type of lensing, the separation between the images from lensing is so 
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small to be resolved by the ground based telescopes. However, the magnification of 
the background star during lensing is observable. In gravitational lensing inside the 
Milky Way galaxy, unlike to the cosmological lensing, the positions of observer, lens 
and source are moving with respect to each other and result is a transient lensing 
effect, so-called gravitation microlensing. In this type of gravitational lensing, the 
lens magnifies the background star in a time scale of one month. 



Fig. 1. A schematic illustration of gravitational lensing with the observer, lens and source. This 
figure is adapted from Refilini 


For the gravitational lensing formalism, we use the simple geometry as shown 
in Figure ©• Here the deflection angle is given by a and it is produced by a single 
lens but we can also extend it to the multiple lenses . Using simple geometrical 
relation, we can write the relation between the position of unlensed source /3 and 
images after lensing 6, as follows: 

6Ds = pOs + aDds, (3) 

where the hat sign represents angle as a two dimensional vector perpendicular to 
our line of sight, Ds is the distance of observer to the source, Dg is the distance of 
observer to the deflector (lens) and Dds is the distance of deflector to the source. 
For the single lens, the deflection angle is given by 

, _ AGM 0 

where it can be written as |d| = 4<I>/c^ and $ is the Newtonian potential. 


( 4 ) 
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Since in the linear theory of General Relativity, the perturbation terms have 
linear property, for the case of having multiple lenses at the same distances from 
the observer, the overall deflection angle can be written as the superposition of 
deflection angles from each lens 

4GMi 0 — §i 

a = y „ „ -;-;—. (5) 

^c^Da\e-e,\^ 

For the case of single lens, substituting equation (jH) in equation ([S]), we can write 
the following equation that relates the position of source to the position of lens: 

02 = /30 + 0 |, ( 6 ) 

where the Einstein angle is define by 

4GM Dds 


q2 


c2 DdDs' 

By solving equation m for single lens, we find two solutions for this equation as 
follows: 




( 8 ) 


In the single lensing, the mapping is radial on the lens plane and the distortion in the 
images results from this radial transformation. The flux of stars on the lens plane 
is approximately constant and it is independent of the position of the lens plane.^^ 
Hence, having larger images from the lensing compare to the unlensed source star 
results in receiving more flux from the source star. The amount of magnification is 
the ratio of the image area to the source area as follows 

I I 

From equation ([8]), we have two solutions for this equation, then the total magnifi¬ 
cation would be the sum of magnifications by two images as 

A = H+ -b A- 

/32 + 261 


A = 




-46%- 


( 10 ) 

( 11 ) 


In this equation the angular distance of source with respect to the lens can be taken 
as a function of time. Assuming a constant relative velocity between the lens and 
source, without imposing any acceleration, the angular location of source in terms of 
impact parameter can be written as the combination of minimum impact parameter 
and relative angular velocity, r as follows 

+ ( 12 ) 

where r is taken a constant value for a simple lens. We can write equation (HB) in 
a simple form by dividing this equation to the Einstein angle, 9e as follows: 

+ 2 


A{t) = 


M\/-b 4 ’ 


(13) 
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where u = /3/6e and equation ([T^ changes to the following form of u^{t) = Mq -I- 
(t —Here Ie is the Einstein crossing time and is defined as tE = Qe /That 
is the typical variation of a microlensing light curve where the source star crosses 
the size of Einstein ring. 

Figure ([2]) shows gravitational microlensing light curve with different impact 
parameters. Smaller impact parameter results in higher magnifications of the light 
curve. 



Fig. 2. Left panel represents the shape of images on the lens plane with point like lenses and 
different impact parameters. The right panel shows the magnification of source stars with different 
impact parameters which is obtained by the overall area of images divided to the area of source 
without magnification. This figure is adapted fronffEU 


The other observable parameter is the angular separation between the images 
and from equation ([81) this quantity is 


| 0 +_ 0 -| = y^ 


m 


je- (14) 

For the case that source, lens and observer are aligned on a straight line, the impact 
parameter is set to zero and angular distance between the images would be in the 
order of 20^. Another useful equation is writing the Einstein angle in simple form 
of 


0B = 


‘^Usch 1 ^ 


Ds 


(15) 


where Rsch is the Schwarzschild radius of lens, x = Dd/Dg- We note that the Ein¬ 
stein angle is a function of lens mass, distance of source and distance of lens from the 
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observer. Using the numerical value of Schwarzschild radius, Rsch = 2.95{M/Mq) 
km, equation (1151) can be written as follows 

es = O.Mlma, (i6) 

Mq lOkpc X 

For the case that lens is located at the middle of distance between the observer and 
source (i.e. x = 1/2) and for a typical lens with = 0.5 Mq, the Einstein angle 
for the the background stars in the Galactic bulge or Large Magellanic Clouds is 
about 


Oe,B ulge — 0.691 TOOS, (17) 

0e,lmc — 0.285 mas, (18) 


which is smaller compare to the resolution of ground based telescopes. 

For calculating the relative angular velocity of source with respect to lens, regard¬ 
ing that all three points of observer, lens and source are moving with respect to each 
other, we choose the position of observer as the centre of coordinate system. The 
relative velocities of the lens-Earth and source-Earth are given by vle = vl — ve 
and vsE = vs — ve, respectively. We set h being a unit vector connecting Earth 
to the source star. A unit vector perpendicular to the line of sight can be made by 
combining velocity and vector n as 

V X n 

|v X h| 


The other unit vector orthonormal to h and 1 is m = n x 1. Substituting 1, the unit 
vector rh is written as 


m = 


V — (v • n)h 

|v X n| 


(19) 


Finally the projection of velocity on plane perpendicular to the line of sight and in 
the direction of rh is given by Vm = v • rh or 


— (h • v)^ 
|v X hi 


( 20 ) 


The relative angular velocity perpendicular to the line of sight is given by ^ = 
Ms ~ Ml as follows 


VSE - (n ■ VSe)^ 
DsIvse X h|2 


(vsE — (vsE ■ n)h) 


^LE - (fi- VLe)^ 
A|vle X h|2 


(vle - (vle ■ n)h). 


( 21 ) 

For the perpendicular component of velocities along the line of sight (i.e. uj.), the 
relative angular velocity simplifies to 


_ vs,j_ — ve.l vl,j. — ve,l 

A • ^ ’ 

Here, we assume that all three points are almost free particles at least within the 
time scale of lensing and velocities are constant vectors. In reality as the Earth is 
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moving around the Sun, observer on the Earth is an accelerating point. We will 
consider this effect, so-called the parallax effect in the microlensing events. Using 
the typical velocity and distance for the source stars in the Galactic Bulge, the 
angular velocity is given by 


= 15.0^as/day 


Ds y 

8.5kpc ) 


vs.j_ — ve,j_-(vl._l — ve,j_) 

X 


220km/s 


(23) 


This velocity depends on the line of sight and vary from the Spiral arms direction to 
the Galactic Bulge direction. Substituting this equation in the definition of Einstein 
crossing time, the time scale of gravitational microlensing events is 


tE = 4:5.6day (—^ f ^ \ ' ( 

\8.5kpcJ \Q.5 MqJ V 


1/2 /i_^xi/2 


-1 


(24) 


In the direction of Galactic bulge the typical time scale of tg is around one month. 
The first microlensing event has been reported in the direction of Large Magellanic 
Glouds by MACHO and EROS collaborations in 1993.1^^11^ Figure ([3]) shows the 
light curve of this event in the red and blue filters. The observation was done in 
two different filters to distinguish the microlensing events from the variable stars, 
noting that gravitational lensing is independent of frequency of light and all the 
light rays follow the same geodesics in the space-time. 

From fitting the observational data of microlensing light curve with the theoreti¬ 
cal light curve in equation (11311 . we can obtain the three parameters of uq (minimum 
impact parameter), tg ( the time of maximum magnification) and ts (the Einstein 
crossing time). The only parameter that has physical information is ts and from 
equation (I24L this parameter is a function of mass of lens, distance of source star 
and lens from the observer and relative transverse velocity of lens and source with 
respect to the observer. In the microlensing experiments the main aim is to iden¬ 
tify the physical parameters of lenses, however the degeneracy between the lens 
parameters makes this goal difficult. However, with the statistical methods for an 
ensemble of microlensing data, we can provide the distribution function for the lens 
parameters. The other methods for breaking degeneracy between the lens param¬ 
eters is with including extra information as the motion of Earth around the Sun 
and finite-size effect of the source stars. We will discus extensively the second order 
effects later. 


4. Statistics of Microlensing events 

The gravitational microlensing studies that initiated by Bohdan Paczynski in 1986 
was for detecting massive compact halo objects of Galactic halo. The idea was 
monitoring huge number of stars in the Large and Small Magenllanic Clouds and 
looking for detection of microlensing events. Finally with counting the microlensing 














March 17, 2015 0:16 WSPC/INSTRUCTION FILE Rahvar-ws-ijmpd- 
March-2015 


10 Sohrab Rahvar 



Fig. 3. The first gravitational microlensing candidate in the direction of Large Magellanic Cloud. 
The light curve has been observer in two blue and red filters. The lowest panel represents the 
ratio of best fit to the microlensing light curve. Since gravitational microlensing is monochromatic 
phenomenon, magnification is independent of frequency of light.I^ 


events we can estimate the contribution of lenses that compose the mass of Galactic 
halo. 

In order to estimate the probability of microlensing events in this direction, 
we can ask what fraction of stars in the Magellanic Clouds undergoes lensing if 
we take a snapshot from the stars in this field . Assuming that a microlensing 
event is effectively happen when a source crosses the Einstein ring of a lens, we 
can assign a cross section for the lenses with the size of ttR'^ similar to the cross 
section calculations in the particle physics, where the Einstein radius is defined by 
Re = DiOe- The observation is performed through a cone to the target stars at 
LMC as shown in Figure (jl]). Assuming density of lenses as a function of distance 
from the observer given by p{z), the number of lenses within a differential element 
of dz is given by 

dN{z) = z'^rtdz, 

^ ' M 

where fl is a solid angle towards target stars and M is the typical mass of lenses. 
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Fig. 4. Observer cone starting from the observer to the source stars at the background (painted 
by the green colour). Lenses along this cone cover the line of sight by the Einstein disks depicted 
by black circles. The ratio of overall angular area covered by the Einstein rings over the overall 
solid angle, represented by Q, provides the probability of microlensing event at a given instant of 
time. 


Multiplying dN{z) to the area covered by the Einstein ring of lenses at the distance 
of z, represents the area of lenses that covers the area of our line of sight towards 
the target stars, 

dA{z) = irRl^^z^ndz. (25) 

By dividing dA{z) to the area of cone at the position of z (i.e. S{z) = flz"^), the 
probability of lensing is given by 

dT{z) =TTR%f^dz. (26) 


We substitute the definition of Einstein radius 


i?i = 


^GMiD, 


c(l - x). 


in equation (l26ll and integrating along the line of sight, the overall optical depth is 


AttGDI 


p{x)x{l — x)dx. 


(27) 


This probability similar to the propagation of particles in a given media in the 
particle physics is called the optical depth. An interesting point about this equation 
is that since area covered by the Einstein ring is proportional to the mass of lens, that 
cancels with the denominator in equation (12611 and the optical depth is independent 
of the mass of lenses and it depends only on the density of lenses in the halo. 
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In order to have an estimation from the optical depth, we calculate the contri¬ 
bution from the spherical structure of Galactic halo, using the density of halo as 
follows 

, , dM 

Substituting this equation in the definition of optical depth in (l?7l) . the optical 
depth is obtained in the order of r ~ (u/c)^. where = GM/Dg is the rotation 
curve of Milky Way driven by the Galactic halo and is in the order of 220 km/s. 
Here we assume that hundred percent of halo is made of MAGHOs. Substituting 
the numerical value for v, the optical depth is r ~ 5 x 10“^ where the target stars 
are in the Magellanic Glouds. This means that observing 10^ stars in this direction, 
almost 5 stars undergoes microlensing magnification. We will see later that from 
results of microlensing experiments, MAGHOs just making less than 20 percent of 
halo mass, then we expect to see only one event out of ten million stars in the 
Magellanic Clouds. 

The other target for gravitational microlensing observations is the Galactic 
Bulge, where this direction has the advantage of large number of bright source 
stars and large column density of stars in the Spiral arms and Galactic Bulge. This 
direction has been used for exploring extra-solar planets orbiting around the lens 
stars. Also one of the important plans in this direction is using microlensing as 
the astrophysical tool for studying the physical properties of the source stars. For 
this direction we can estimate the optical depth similar to the direction of Galactic 
halo. Here stars in the disk of Galaxy play both the role of lens and source stars. 
Assuming that disk has cylindrical symmetry, the mass of disk as a function of 
distance from the centre is given as p{r) = dM/{2Trhrdr) where h is the width of 
the disk and r is the distance from the centre of Galaxy. Substituting this density 
in equation (l27l) . the optical depth obtain about 

2GMdisk 

^ “ c2/i ■ 

For the case of stellar mass with Mdisk h = 300 pc, we get the value of 

T ^ 3 X 10“®. A realistic value for the optical depth needs detailed calculations in 
different directions of Milky Way galaxy .1^1^ 

The other statistical parameter for counting the microlensing events is the rate 
of events per a given time (traditionally given in terms of year). Let us assume that 
lenses are moving with constant transverse velocities and we monitor a patch of 
sky for the duration of Tabs- In this case, each lens spans a rectangular area at the 
position of z form the observer. The total area spans by lenses is 

dA{z) = 2REvTobs X ^nz'^dz. (28) 

Ml 

We multiply the denominator and numerator of equation (E51) to Re, moreover 
multiplying this equation to the number of background stars. The result is the 
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number of microlensing events for a given observation time of Fobs- Then, the rate 
of events per number of background stars and observation time is given by 



(29) 


Comparing this equation with the definition of the optical depth in equation (I26L 
the rate of events can be written in the simple form of 



(30) 


or 



(31) 


TT tE 


This equation provides the rate of microlensing events in an ideal condition where all 
the ongoing microlensing events are observable for a telescope. In practice we miss 
very short duration microlensing events as well as events longer than the duration 
of experiment. The efficiency of detection of an event with the duration of is 
given by e(t£;) function and provides fraction of events that can be detected by a 
telescope. This efficiency function obtain by synthetic observation, using a Monte- 
Carlo simulation and producing catalog from the microlensing data in the computer. 
Now, let us assume an experiment that detects N microlensing events. Each event 
has an Einstein crossing time of where i represents the index of event. The 
corrected rate of events in terms of efficiency function is given by the inverse of 
this function and divided to the exposure time and Number of background stars as 
follows 



(32) 


r = 


The other important parameter that can be compared directly with the theory 
is the optical depth. We can express the geometrical definition of optical depth in 
terms of temporal definition, by means that optical depth is equal to the overall 
duration of microlensing events divided to the exposure-time times the number 
of background stars. The overall Einstein crossing time has to be corrected by the 
efficiency function, by means / e(t^E ) ■ Since the duration of events depends 

on the impact parameter, for the case of zero impact parameter the duration is 2tE 
while for the impact parameter with mq = 1, according to the definition of an event 
with Amax > 1.34, the duration is zero. For the impact parameter ranging form 
zero to one, the average impact parameter is 7r/2 x tE- Then, the temporal optical 
depth is defined as 



(33) 


In what follows, we will compare the observed optical depth with the theoretical 
value. 
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4.1. Observational results from detection of MACHOs 

The microlensing experiments of EROS,l^l^ MACHO, MOAI^ and OGLEP 
have started observation of Large and Small Magellanic Clouds in 1990s for search¬ 
ing microlensing signals from compact objects as candidate for dark matter in the 
halo of Milky Way galaxy. These telescopes were observing in survey mode by moni¬ 
toring a large area of target stars and investigating the light curves for microlensing 
events. There were also small size follow-up telescopes for making full coverage of 
light curves with high cadence and higher signal to noise photometry data points. 
This mode of observation is still in used for extrasolar planet observations with 
microlensing method. 

After producing photometric images, the main part of observation is produc¬ 
ing light curve from the images. Especially in the crowded fields it is difficult to 
distinguish the time variation of background stars from the microlensing events. 
The observation groups developed image difference method where change in the 
flux of a star in the crowded fields can be identified by differencing images from a 
template image of the field.l^ The result is identification of time variation of a star 
as a function of time. Due to variability of source stars, there are large number of 
backgrounds that has to be cleaned from the signals. In order to select microlensing 
signals, stable stars are chosen for flux variation analyzing and in the direction of 
Magellenic Clouds, out of millions of background stars few microlensing candidate 
has been found by the observation groups. 

The next step is analyzing light curves by fitting simple microlensing light curve 
and extracting the parameters of events. The only physical parameter in the light 
curve of single lens is the Einstein crossing time which depends on three parameters 
of lens mass, distance of lens and source from the observer and transverse velocity of 
lens. Due to degeneracy in the parameters, we need to perform statistical analysis in 
terms of number and duration of events. Since the number of events in the direction 
of LMC/SMC is small, there would be large Poisson error. The final number that 
will help us to estimate the contribution of lenses in the Galactic halo is from the 
optical depth calculation. An extensive discussion on the optical depth calculation 
from the experiment can be found in Ref.^^ The most challenging problem in the 
optical depth calculation is the blending effect. This effect results from the mixture 
of background stars within the Point Spread Function (PSF) of source star. 

The blending effect can have two important effects of (i) modifying the shape 
of light curve that provides a wrong estimation for the value of Einstein crossing 
time of a microlensing event and (ii) we get wrong estimation for the number of 
background stars in the field. In the blending effect out of many stars inside PSF 
one of the stars in PSF disk is magnified, then the overall flux is given by 

F{t) = FoA{t) + F', (34) 

where Fq is the microlensed star and F' is the overall flux of background stars. 
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Then the observed magnification is given by 

FoA{t) + F' 


Aobsii) — 


Fo + F' 


We can rewrite this equation in the following form of 

Aobsit) = bA{t) -f (1 — 6), 


(35) 


(36) 


where b is the blending parameter and is given by 6 = Fq/{F' + Fq) and ranges be¬ 
tween zero for when there is maximum blending and one when there is no blending. 
For the case of 6 < 1, the observed magnification is always smaller than intrinsic 
magnification (i.e. Aobs < A). This means that for a given impact parameter of 
ito, the blended microlensing event will suppressed and the result would be smaller 
Einstein crossing time for a microlensing event. One the other hand since oc y/M, 
then estimated mass for the lens would be smaller than real value. 

The other problem of blending is that we will get an incorrect number for the 
source stars in the background. Since the optical depth from the experiment in 
equation (1331) is given in terms of the number of source stars, underestimating back¬ 
ground stars will provide overestimating the numerical value for the optical depth. 
The microlensing experiments updated their measured optical depth by discarding 
some of their microlensing candidates as well as having better value for the number 
of background stars from blending estimation. The updated value for the optical 
depth of MACHO group in the direction of LMdSlI^ is r = (0.4 ±0.1) x 10“^ 
after 5.7 years of collecting datcP^ while the first estimation was a larger number of 
T = 1.2)^q’3 X 10“^ . The EROS collaboration also from 6.7 years data collectiorP^ 
found three microlensing candidate towards LMC and obtained the optical depth of 
T = (0.15 ± 0.12) X 10”^. In order to decrease the blending effect, EROS group also 
used bright stars for estimating the optical depth. In the direction of LMC, for the 
bright sources there was one microlensing candidate and using this subsample they 
obtained the following limit of r < 0.36 x 10“^ within 95% of level of confidence. 
Figure ([5]) represents an upper limit to the contribution of MACHOs in the Galac¬ 
tic halo. Using the Standard model of Galactic halo,l3^ the expected distribution of 
events in terms of mass of lenses is in the upper panel and excluding halo made of 
MACHOs for a given lens mass is at the lower panel of Figure ([5]) This analysis 
shows that within the range of 10“® to few solar mass MACHOs are making less 
than 20% of halo mass. In this analysis the mass function of lenses are taken as 
Dirac-Delta mass function, given in the X-axis of Figure ([5]). 

The other method for estimating blending parameter is using images from the 
Hubble Space Telescope (HST) to compare some parts of sub-fields of target stars 
with diffraction limit images The other possibility is using Lucky imaging 

method to remove the turbulent effects of atmosphere and producing least blended 
images from the field of stars.^^^^ 
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Fig. 5. For the Standard model of Galactic halo, EH the expected distribution of events in terms 
of mass of lenses in the upper panel and excluding halo made of MACHOs for a given lens mass 
is at the lower panel. This Figure is adapted from Ref.E^ 


5. Second order effects in single microlensing events 

We have seen that for the case of microlensing events with the single lens, the 
Einstein crossing time which is the only physical parameter of lensing, depends on 
the mass, distance and transverse velocity of lens with respect to the line of sight. 
Due to this degeneracy we can not identify the parameters of an individual event. 
However, as we will see in this section, the second order effects as parallax and finite 
size of source star will help us to break partially the degeneracy between the lens 
parameters. 
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5.1. Parallax Effect 

In simple microlensing light curve, the relative motion of observer-lens-source in 
equation (1221) is non-accelerating. Any accelerating motion of these three points 
results in deviation of light curve from simple microlensing light curve. The simple 
microlensing light curve is symmetric as shown in Figure ([2]). The annual motion 
of Earth around the Sun or rotation of lens around companion or source around 
companion results in deviation of light curve from single microlensing light curve. 
Figure ([6|) demonstrate the rotation of Earth around the Sun and subsequently the 
relative motion of observer-lens-source.^ 



Fig. 6. Relative motion of observer with respect to the lens on the lens plane. In the upper panel 
the line of sight to the source star crosses the lens plane while Earth is moving around the Sun and 
this intercepted line follows an ellipse on the lens plane. On the other hand lens follows a straight 
line. Here the position of Sun is located at the centre of ellipse and relative distance between the 
line of sight to the source and lens enters in the single magnification formula of equation dlsll. 
Figure is adapted from Ref.[ini 


In order to calculate the relative velocity, we first calculate the coordinate of all 
components of observer, lens and source, then derive velocity from the equation of 
motion. The motion of earth around sun, projected along the line of sight is give by 

fE(t) = a© cos^(t)i-I-a© cos/3sin^(t)j, (37) 

where Sun is located at the centre of coordinate and a© is the orbital radius of 
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Earth, ^ {t) = cot+ (j)o is the orbital phase of Earth, uj is the orbital angular velocity 
of Earth and /3 is the angle between the ecliptic and deflector plane. Here we assume 
negligible eccentricity for the earth orbit. Hence the velocity of earth from derivation 
of equation (l37l) is given by 


VE{t) = —ia^uj + ja^^uj cos f3 cos ^(t). 

The position of lens projected on the lens plane can be given as follows: 

= iRs 


t-to 
tE 


cos 9 l + uo sin0 


■ jRe 


t-tp 

tE 


sin 9e + Uo cos 9 


(38) 


(39) 


where 9l is the angle of trajectory of lens with respect to the semi-major axes of 
earth orbit, projected on the lens plane and 9l G [0, 27r]. We note that in order to 
cover all the possible configurations, the velocity of lens needs to be both positive 
and negative signs. Hence we let tE has both positive and negative signs. The time 
derivation from equation (IMl) results in velocity of lens as follows 

VL{t) =i^cos9L+j^sm9L, (40) 

tE tE 

where vl = RE/tE is the transverse velocity of lens (and can has positive and 
negative signs). Similar argument can be used for the motion of source star. The 
velocity of source is given by 


us{t) = ivs cos9s + jvs sin0s, 


(41) 


where vs is the projected velocity of source star parallel to the lens plane and 9s is 
the angle with respect to the semi-major axis of Earth orbit. The apparent relative 
velocity of these three points is given by equation (1^ . In more complicated case, 
in addition to the observer, the source star and lens can undergo an accelerating 
motion, if they belong to gravitationally many-body bounded system. The case for a 
binary lenses and binary sources can be found in Ref.^^^ For simplicity, we assume 
that the motion of lens and source are non-accelerating and contribution from the 
annual motion of earth introduces time variation in the relative velocity of these 
three points. 

In order to calculate the magnification of source star from equation (m , we 
need the relative distance between the observer-source intercepting line in the lens 
plane with the position of lens. We use equations dSIl) and dSi) to calculate this 
distance. In order to calculation magnification from equation (11311 . all the scales 
should be projected on the lens plane and normalize to the Einstein radius. The 
relative distance on the lens plane is given by u(t) vector as follows 

u(t) = i ((5rtcos^(t) — -—— cos 9 l — uosin^i 

V tE 

+ j ^(5MC0S,5sin^(t) — ^ ^ sin0i — uo cos9l^ , (42) 
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where Su = a©{l — x)/Re- In literature, 5u is also expressed by the symbol of 
also the projected velocity of lens on the observer position, v = u/(l — x). These 
two parameters are related to each other by 


V = 


Qe 

tsSu 


(43) 
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Fig. 7. First microlensing event with parallax signature. The effect is deviation of light curve from 
symmetric simple microlensing shape. The upper panel is magnification verses time in R-band and 
lower panel is at B-band. The parameters of lens are tE = 110 days and projected velocity of lens 
on solar position u = t?^/(l — x) = 75 zb 5 km/s or Su = 0.21. Figure is adapted from Ref.^^ 


Figure © shows the first parallax microlensing event observed by MACHO 
group The parallax effect causes simple microlensing light curve deviates from its 
symmetric shape. For the case of Figure ©, the Einstein crossing time is ~ 110 
days, comparable to one year annual motion of earth. Having longer events, the 
trajectory for the impact parameter, u(t) follows more curved path compare to the 
short duration events. Hence for events with <_e in the order of hundred days, the 
parallax effect is more effective. 

From equation (|42|) . the extra parameters that enter in calculation of u(t) are 
Su and 0l which can be obtained from fitting theoretical light curve (including the 
parallax effect) to the observed data. From the five parameters of uq, to, Ie, 0l and 
Su enters to the magnification function, only Ie and Su have physical information. 
From simple lensing, we have seen that knowing Ie, we obtain a constrain between 
the mass, distance and relative velocity of the lens from the observer. The parallax 
parameter, Su has an extra information and from equation (1421) . Su can be written 
as 


Su = 


1^Z?^)-1A(L^)1/2^ 


(44) 
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Fig. 8. Relation between the mass of lens (in solar mass unit) verses distance of lens from the 
observer in terms of a; = Di/Dg. The parameter of Su = 0.2 is assumed to be measured from the 
microlensing light curve. Here we set source star is located at Da = 8 kpc in the Galactic bulge. 


unlike to the Einstein crossing time this parameter depends only on the mass and 
distance of lens from the observer. Figure ([8]) shows constraining relation between 
the mass of lens and relative distance of lens and source for a given value oi 6u = 0.2. 
In Figure dl]), we assume source is located at Galactic bulge, where ^ 8 kpc. We 
can estimate the mass and distance of lens from the observer if we use models for 
the distribution of mass in the Milky Way galaxy and the mass function of lenses.!^ 


5.2. Simultaneous parallax measurement by a ground-based 
telescope and a telescope orbiting around the Sun 

Traditionally, the displacement of angular position of an object compare to the 
reference points at infinity by two different observers is also called parallax effect. 
This method has been used in the astronomy for measuring the distance of remote 
stars by measuring the displacement of the angular position of a given star relative 
to the background stars at two different times, preferably with six month difference. 
We use this technique for observation of gravitational microlensing events, but since 
microlensing is a transient event, we need to perform observation at the same time 
while from two different positions. 

The degeneracy between the lens parameters in the microlensing observation 
has been known since 50 years and for the first time Refsdal and Liebes proposed 
breaking this degeneracy by observation of microlensing events from two different 
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positions. One of the observers is located on the Earth and the other one is a 
satellite telescope, orbiting around the Sun.I^^JEH Let us assume that the projected 
distance between the Earth and satellite on a plane parallel to the lens plane, at the 
position of solar system is d±. From the observation of light curve we can identify the 
minimum impact parameters as well as the time difference between the maximum 
peak of the light curves for the two observers. Lets us use ue for minimum impact 
parameter of the Earth observer, us for minimum impact parameter for the satellite 
observer and At for time difference between the maximum magnification for the two 
observers. The projection of d± on the lens plane normalized to the Einstein radius 
is related to the difference in the impact parameters and time difference between 
the peaks as follows: 


iSUf = ( — f + {uE-usf, (45) 

tE 

where 51/ = (1 — x)d±/RE. The geometrical details are depicted in Figure (|9]). The 
relation between SU and 6u or tte that is introduced in previous section is given by 
SU = tte X- d±/a(^. Hence the parallax parameter in terms of two components in 
equation (H51) is given by 


a© At 

Ti'E = 3 —( —, me — us), 
a± lE 


(46) 


where At/tE is along x-axis, parallel to the trajectory of the lens and y-axis is 
perpendicular to it. 

From Figure (jH]), we note that since u ~ t/tE (i.e. impact parameter is propor¬ 
tional to time over the typical leasing time), then |ue — us| — At/tE and we would 
expect that two components of tte vector to be in the same order.lSSl Also we can 
identify the trajectory of source star in the lens plane by tana = [ue — us)tE/At. 
Recently it has been possible to observer simultaneously the microlensing events 
from the Earth and Spitzer Satellite orbiting around the Sun.^Sl Dedicated tele¬ 
scopes orbiting around the Sun in future will break partially the degeneracy of 
microlensing parameters. Using tE form the light curve of event and the average 
value for the relative transverse speed of lens, we can estimate the Einstein radius. 
Combining with equation (|4^ , we can calculate the mass and distance of lens from 
the observer. The main impact of this kind of projects would be identifying the dis¬ 
tance of lenses, their masses and more impotently the distance of planets orbiting 
around the lens stars. 


5.3. Terrestrial Parallax effect 

The other kind of microlensing parallax observation is the terrestrial parallax. The 
parallax observation from Earth has been used for long time for measuring distance 
of close-by objects such as moon. However, the accuracy of parallax observation 
from the earth for distant objects is very low and practically it is useless. In the 
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Fig. 9. Relative distance of Earth and a satellite orbiting around the sun from the position of 
lens. The ellipse is the projected trajectory of earth on the lens plane, and ug are the minimum 
distance of earth and satellite with respect to the lens, respectively. 5U is the distance of the Earth 
and satellite normalized to the Einstein radius at the maximum magnification for light curve of 
each observer. 

case of very high magnification microlensing events, observers at different positions 
on the Earth can see the peak of light curve at different times. Using this argument 
that Am ~ At/tE, the time difference between the observation of peaks for two 
observer would be At ~ x Am. Taking ~ 10 days, Su = D^/Re where 

~ 10^ km is the diameter of earth and Re ~ 1 A.U., we get At ~ 1 min. The 
first terrestrial parallax microlensing event has been observed in OGLE-2007-BLG- 
224 event where the reported time delay in peak of microlensing light curve between 
three observatories of South Africa, Ganaries and Ghile were about one minute.l^ 

5.4. Finite source effect 

The other important perturbative effect compare to the simple microlensing light 
curve is the finite-source effect. The finite-source effect happens when the size of 
source projected on the lens plane is not negligible compare to the impact parameter 
of lens from the centre of source star. Then the different parts of source will magnify 
with different factor and result would be deviation of light curve from the simple 
lensing. 

Let us assume that flux of star is also non-uniform while disk of star is rotation- 
ally symmetric around the centre of disk. We choose the coordinate at the centre of 
disk of source star, using the magnification factor depends on the distance of lens 
from each point of source disk. The configuration is shown in upper panel of Figure 
(HUD and the overall magnification factor is given by 

^ ^ iZo Ivlo A{w)F{v)vdvde 

Ie=o Ivlo Fiv)vdvde 


( 47 ) 
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where w = (u^ + v'^ + 2uv 00^9)^/"^ and p* is the size of source star projected on 
the lens plane, normalized to the Einstein radius as shown in Figure (1101) . Detailed 
analytical integration of equation (1471) for a source with constant flux over its sur¬ 
face, is calculated by Witt and Mad^ in 1994. In this case we can simplify two 
dimensional integral to a one dimensional integral. Let us change coordinate system 
from (u, 6) to {w, (f) where (j) is angle between w and a line connecting lens to the 
centre of source star and is given by 

m = arccos)-). 

2uw 

Then the differential element for area is given by wdwdcj). In this case for three cases 
oi w > u, w < u and w = u, the magnification factor is given by 


Mu) = —y 


rl“+P*l 


Trp: 

TT 


'* J\u-p^\ 


V4 + 


arccos(---) aw 




2uw 


+ - u)(,p* - m)v^(p*-m)^+4, 


(48) 


where H is Heaviside function and is equal to one if u < p*, otherwise it is zero. The 
maximum magnification in finite-size effect at u = 0 is given by A^ax = -^/l -b 4/pJ 
which is larger for the small size stars. For the case of p* —> 0, we will have singularity 
in the maximum magnification. 

Figure (flUl) (lower panel) shows a sample of light curve with finite source effect 
for a very high magnification microlensing event .1^ One of the advantages of high 
magnification events is that lens can probe the luminosity of surface of source star. 
Since the luminosity due to limb darkening is not uniform and the apparent lumi¬ 
nosity of stars decreases around the edge of star, the gravitational microlensing is a 
unique tool to probe the limb-darkening of stars few thousand parsec far from us. 
A simple limb-darkening model is / = /q [1 — wa( 1 ~ cos<())], where (j) is the angle 
between the line of sight and a vector normal to the surface of the star. The best 
value of limb-darkening parameter in Figure (flOl) for event OGLE-2004-BLG-254 is 
un = 0.70 ± 0.05 and uj = 0.55 ± 0.05 where the size of star is p* = 0.0400 ± 0.0002. 
The parameters of limb darkening can be measured for any high magnification event 
with a perfect coverage of the peak of light curve. One of these samples is OGLE 
2008-BLG-29(P 


5.5. Detection of stellar spots in Finite size effect 

The capability of gravitational microlensing as a natural telescope can also be used 
for probing the stellar spots on the surface of stars located at the other corner of 
Milky Way galaxy.l^l^ The presence of spot produces perturbations in the mi¬ 
crolensing light curve when lens get closer to the spot. We expect to observer such 
an effect in very high magnification events when the impact parameter is in the 
order of projected size of star on the lens plane. In this case the magnification on 
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HJD-2450000 


Fig. 10. Upper panel demonstrates the distance of projection of source star on the lens plane 
from lens. The lower panel is a sample of very high magnification event OGLE-2004-BLG-254 and 
the finite source effect on the light curveE^ 
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time [source-radius crossing time 


Fig. 11. Simulation of a microlensing light curve with zero impact parameter and a spot with the 
size Vs = 0.2, normalized to the radius of the star. The solid line and dashed line are with different 
distance of spot from the lens and dotted line is the microlensing light curve without spot. This 
figure is adapted from RefP! 


the source star is given by 

/g F{r')A{\r - r'|)dS' - ^ Fn{r')A{\r - r'|)dS' 


(49) 


where r is the distance of centre of source from the lens and r' is the position of 
points on the surface of source from the centre of source disk. F{r) is the spotless 
flux of source and Foir) = F{r) — Fs{r) is the flux increment between the spotless 
source and flux of spots. The integral index of S represents overall area of source 
and the index of S' represents integration over the spot. In the case that lens crosses 
the position of spot on the lens plane we can estimate the time scale of perturbation 
by t, ~ / X p* X where / is the ratio of spot size to the source size. In the case 
of binary lensing where the singular area of lens unlike to the single lens is a curved 
path, the probability of caustic crossing of spot is higher than the case of single 
lensing .1^ 


5.6. Self-Lensing and observation of compact companions around 
source star 

Almost two-thirds of stars are in binary systems. For the case that these systems 
locate at the edge-on position with respect to an observer, the companion star can 
block light receiving from the other star during eclipsing and result is decreasing the 
overall light from the two unresolved stars. This method has been used for long time 
to identify binary systems. Using the spectroscopic as well as photometric data, one 








March 17, 2015 0:16 WSPC/INSTRUCTION FILE Rahvar-ws-ijmpd- 
March-2015 


26 Sohrab Rahvar 

can find the parameters of edge-on binary systems. 
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Fig. 12. Light curve of self-lensing form Kepler Photometry KOI-3278. The left panel is the 
microlensing signal during the self-lensing and right panel is the occultation of white dwarf. The 
residual between the model and data also is written in the lower panel. The binary system is made 
of sun-like star and a white dwarf. This figure is adapted form Refpa 


After evolving one of companions in a binary system, it can turn to a black hole, 
neutron star or white dwarf. In this case, instead of blocking light of companion star 
during eclipsing, the compact companion star can magnify light receiving from the 
other star. Then, we would expect having a small magnification due to self-lensing 
in the binary system. The first prediction of this effect has been done by Andre 
Maeder in 1973, long time before the first detection of microlensing event.For the 
short period binary systems, we would expect to observing repeating self-lensing of 
companion star by its companion. With developing ground based and space based 
telescopes, it was predicted to observe this phenomenon in practice.^^^^ 

The first case of self-lensing event in binary system has been discovered in 
database of Kepler space-telescope in 2014.1^ This event, so-called KOT3278 is 
a binary system composed of sun-like star and a white dwarf orbiting it every 88.18 
days. The microlensing transit occurs in 5 hours and the result of lensing is 0.1% 
magnification in the flux of companion star. Figure m shows the first self-lensing 
candidate of KOT3278 in Kepler data. 


6. Astrometric and Polarimetric Observation 

In this section we will introduce two alternative methods of astrometric and polari¬ 
metric observations of microlensing events. In the astrometric observation, we are 
interested in to detect the position of images or in another word the displacement 
of source position due to dynamics of images during the lensing. The other method 
is the polarimetric observations of microlensing events for those source that have 
polarized signal. These two methods will produce extra-informations and enable us 
to break completely the degeneracy between the lens parameter. Moreover, from 
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the polarimetric observation we can extract extra information about the nature of 
source stars. The observation of these two parameters will open a new window to 
study the physics of stellar atmosphere and application of microlensing in the stellar 
astrophysics. 


6.1. Astrometry of microlensing events 


In the previous sections we have discussed about the photometry of gravitational 
microlensing events where the flux of light receiving from a source star becomes 
brighter during the transit of source over the Einstein ring. The only physical pa¬ 
rameter that we can extract from the photometric observation is the Einstein cross¬ 
ing time which depends on the mass, distance and velocity of a lens. In this section 
we discuss on possibility of astrometric observation of microlensing events. 

During microlensing by a single lens, two distinctive images produce at the either 
sides of the lens position. The separation between the images from equation (H and 
m is in the order of Einstein angle and it is about ~ 1 mas. This angular distance 
between the images is smaller than the resolution of ground based telescopes. High 
resolution space-based telescopes or adaptive optics may resolve these images in 
future.l^i^^ In the astrometry we can measurement of the displacement of centre 
of light during the lensing which do not need extremely high resolution to resolve 
milli-arc second separation between the images from the lensing. Here, we need to 
identify and follow the centre of Point Spread Function (PFS) of the images during 
the lensing. The centre of light of images with respect to the position of lens is given 

by 


\A+\e+ + \A-\e- 
l^+l + |A-| 


(50) 


where |H+| and \A~\ as discussed before are the absolute magnification of two 
images at the angular positions of 0+ and 9~. Substituting the magnification and 
position of images from equation ([5]) results in 


0c= P 


30| 


/32 


20 |' 


(51) 


The observable parameter in the gravitational microlensing is the displacement of 
the centre of light of source star during the lensing. This parameter can be obtain 
hy 60 — Oc — P and the result is 


59 = P- 


92 


20|- 


(52) 


Figure (fT^ represents the absolute value of astrometric shift during microlensing 
for an event with the minimum impact parameter of Po = 0.59e- The maximum 
displacement of centre of light from d\69\/d\P\ = 0 happens at /3 = V^9e which 
results in the maximum displacement of centre of light, SOmax = The 

observation of this quantity directly provides the observation of Einstein angle. 
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We note that while the maximum photometric magnification happens at the 
smallest impact parameter, the maximum astrometric shift happens before the 
source star enters the Einstein ring. The other property of astrometric shift is that at 
the large impact parameters, 69 falls as 1//1 while for the magnification 5A oc l/zS"*, 
means that we can detect the astrometric shift even at the larger impact param¬ 
eters. This property results in a larger astrometric optical depth compare to the 





Fig. 13. Astrometric shift of images in microlensing event with the minimum impact parameter 
of /3o = 0.5. The impact parameter and astrometric shift angles are normalized to the Einstein 
angle. 


photometric optical depth.^^^^ For the case of GAIA satellite where the accuracy 
in the astrometric shift is in the order of micro-arc second, the optical depth can be 
obtain with the accuracy limit of telescope. In this case from (l5^ . assuming 5t as 
threshold for the observation, the corresponding maximum impact parameter ob¬ 
tain as /3 = Q\I6t- Comparing the astrometric optical depth with the photometric 
optical depttP^ shows that in the former case optical depth is about 10“® while for 
the later one, it is in the order of 10“^. 

In the case of accompanying the astrometric observation with the photometric 
observations, having high cadence and high photometric precision, it would be pos¬ 
sible to detect the perturbation effects such as parallax or finite source size effects. 
This will help us to break completely the degeneracy between the lens parameters.!^ 
Let us assume accompanying astrometric observation with the parallax measure¬ 
ments in the light curve from the ground based of space based observations. In this 
case we obtain the parameter of 6u which is given in terms of Einstein angle and 
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distance of lens and source as follows 


Su = 


Qe(l - x) 
Dsx9e 


For the microlensing events in the direction of Galactic Bulge the distance to the 
source stars is almost 8 kpc. On the other hand from astrometric observation, we can 
measure 9e- Substituting these parameters in the parallax equation, we can extract 
X parameter which is the relative distance of lens to the source. On the other hand, 
using the equation for the Einstein angle and distance of lens and source, we can 
obtain the mass of lens and subsequently the transverse velocity of lens-source from 
measuring the Einstein crossing time. 

The other helpful observation is the finite size effect of source star which enables 
us to measure the projected radius of source star on the lens plane, normalized to 
the Einstein radius (i.e. p* = xRg/R e)- The size of source star, after implementing 
the extinction correction of interstellar medium, can be obtained from the position 
of star in the colour-magnitude diagram.^ One of the problems with identifying 
the position of the star in the colour-magnitude diagram is the blending effect of 
either lens star or stars in the background that can alter the position of star. 

Measuring p* and Rs can put constrain on p* = RsX/Re and the result is 
constraining between the mass of lens and distance from the observer. On the other 
hand from the astrometry we can measure the Einstein angle 9e which combining 
with the finite size effect, we can measure the distance of source star from the 
observer. 


6.2. Polarimetry observation of microlensing events 

The scattering of photons by the atmosphere of stars produces linear polarization. 
The mechanism of producing polarized light is due to scattering of light by free elec¬ 
trons of atmosphere in the hot stars,^ scattering of light by atoms and molecules 
in the main sequence starJ^ and scattering by the dust grains in the cold stars.® 
The effect of polarization in gravitational microlensing has been studied in Super¬ 
nova explosions where local polarization on the gas of a symmetric explosion can be 
magnified by another intervening compact object as a gravitational lens.® The ap¬ 
plication of this phenomenon to the microlensing events within Milky Way galaxj® 
is aiming to enhance polarization signals from the source stars during lensing. This 
method enables observers to study the atmosphere of stars thousands of parsec far 
from Earth, somewhere at the Galactic centre or Spiral arms of Galaxy. 

Assume a symmetric polarization of source star due to the scattering of light 
by the atmosphere of a star, in the gravitational microlensing, different parts of 
source disk magnify with different amount and for the smaller impact parameters, 
the higher is the magnification. The result is breaking symmetry of polarization 
pattern and producing a net polarization signal from the source star. The extra 
information from the amount of polarization and time dependence of it during 
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lensing can break the degeneracy between the lens parameters as well as provide a 
unique information about the atmosphere of source starsP^^ 



Fig. 14. Polarization ellipse with semi-major axis of A and semi-minor axis of B. The semi-major 
axis is rotated with the angle of 6 in counterclockwise direction with respect to the coordinate 
system. This Figure is adapted from referenceE^ 



Fig. 15. Projection of source star on the lens plane where Z direction represents the line of sight 
from the observer to the source. The intensity of lens is decomposed into radio direction (i.e. Ii) 
and transverse direction (i.e. Ir)- Here we adapt the convention by Chandrasekhar.^^This Figure 
is adapted from Reference.!^ 
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The polarization of a source star is given in terms of four Stokes parameters 
which is expressed by Sj , S'q, Sy and Sy- Let us assume electromagnetic filed in 
the directions of X and Y axis perpendicular to the propagation of light, varies as 
E'x = klexp(ia;t) and E'y = B exp{iujt + 7r/2). Now we rotate coordinate system 
along the propagation of light, in clockwise direction by the angle of 9 (as shown 
in Figure ITdl) and in the new coordinate system calculate the Stokes parameters. 
While the circular polarization for stars with taking into account the stellar spot 
is non zero,l^^ for simplicity we set Sy = 0 for stars without spots. The Stokes 
parameters on the surface of star with the area of S defined as follows: 


Si= [ <\Ex{s)\^ + \EY{s)\^>d^S 

JS 

= [ I{9)dPS, 

JS 


{A{s)^ + B^{s)) (fS 


(53) 


where /(/i) represents the intensity of star in polar coordinate on source disk and 
/r is the cosine between the line of sight and vector normal to the surface of star. 
Using the Chandrasekhar convention, for the polar coordinate on the surface of a 
star, the intensity of light in the radial and tangent directions are taken as Ii = 
and Ir = B^, see Figure (IISl). 

The other Stokes parameters are 


= ^ < Ej,{s) - EUs) >d?S = J {Aisf 

= —J I-ifJ.) cos29d^S, 


Su = - f < E*x{s)Ey + E^{s)Ex{s) > d^S 
JS 

= J I-in) sm26d^S, 


B^{s))cos2ed^S 


(54) 

— J [A{s)^ — B'^{s)) sin20d^«S' 

(55) 


where I- = B^ — A^ = Ii — I^- 

The apparent intensity of stars for an observer can change due to limb darken- 

in^H 


/(M) = /o(l-Ci(l-Ai)), 


(56) 


where Iq and ci are the parameters for limb darkening function and fi as described 
before represents the distance from the centre of source disk (e.g. /r = 1 is the centre 
and fi = 0 represents limb of star). The high magnification microlensing events can 
probe the parameters of limb darkening.!^ The polarized component of source stars 
also can be given by 

/_(/r)=/oC 2 (l-/r), (57) 

where C 2 represents local polarization fraction on surface of star. The parameters of 
limb darkening and polarization depends on the physics of atmosphere for a given 
star.l^ For the hot stars, using the light scatteting by the free electrons in the 
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atmosphere, Chandrasekhar (1960) calculate the limb darkening and polarization 
from the surface of this type of stars as follow^^ 






1 -f 16.035/r-f 25.503At2 

1-f 12.561^-f 0.331^2 ’ 

0.1171-f 3.3207^-f 6.1522/^2 
1-f 31.416Ai-h74.0112/i2 ~ 


(58) 


For the case of main sequence stars such as Sun, the polarization is mainly due 
to the Rayleigh scattering by neutral hydrogen and partially due to Thomson scat¬ 
tering by free electrons. The polarization of light for the case of Sun is a wavelength 
dependent function and it can be modelled by using the structure of atmosphere of 
Sun. The best fit to the theoretical model is given 


where parameters with the index of A depends on the wavelength. For the case of 
cool stars, depending on the atmospheric model and ingredients of dusts, the polar¬ 
ization function can change.!^ Using the magnification factor in the gravitational 
microlensing, the Stocks parameters can be written as follows: 


/Sq 

Su 

\Si 


= p* 


p dp d9A(u) 


' —I-{p,) cos20' 
I- (/i) sin 20 

lip) 


(60) 


where u is the impact factor (distance of lens from each element on the source star) 
and A{u) is the local magnification factor. 

It is supposed that the next generation of microlensing experiments will perform 
polarimetry as well as astrometry observations beside the photometry observations. 
This supplementary observations will open a new window for scanning the atmo¬ 
sphere of remote stars and investigate the details of their structures. Moreover, 
using this method, we can study defects on the surface of stars as stellar spots and 
even measure the strength of their magnetic field. 


7. Summary 

In this article we reviewed the astrophysical applications of single-lens gravitational 
microlensing events. The first microlensing candidate has been observed almost two 
decades ago and the aim of microlensing experiments at that time was discovering 
the Massive Astrophysical Compact Halo Objects (MACHOs) in the galactic halo. 
However, the suggestion of exploring planets in the binary lenses with a parent star 
and planet orbiting around it, opened a new window to the microlensing observa¬ 
tions and since more than ten years, microlensing experiments are mainly focused 
on exoplanet detections. With this method more than ten exoplanet at the distance 
of snow-line from the parent stars have been discovered. Moreover, using large field 
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telescopes and follow-up observations, more than 3000 microlensing event have been 
discovered each year in the direction of Galactic Bulge. 

Using dedicated space based telescope by means of observing events both from 
Earth and space, it would be possible to do both parallax measurement and astro- 
metric observations. This kind of next generation microlensing experiments could 
measure the mass of lenses, as well as distance of lenses from the observer. Since 
part of lenses are too faint to be observed or they are compact objects as black hole, 
neutron stars or white dwarfs, it would be possible to produce the mass function of 
stars and stellar remnants. Also knowing the distance of stars, it is possible to study 
the structure of the Milky Way and distribution of stars along the line of sight. 

The other phase of microlensing observation would be the astrophysical appli¬ 
cation of this method as a natural telescope to scan the surface of source stars. 
This application can be done by photometric and polarimetric channels. With this 
method one can detect small defects as the stellar spots on the surface of stars even 
located at the Bulge of Milky galaxy. Regarding vast applications of microlensing, 
it seems that this technique will be an important astrophysical tool to study the 
stellar physics as well as the structure of Milky Way galaxy. 
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